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1. B-physics and lattice QCD 

For the plenty of beautiful results from recent 
and current B-physics experiments [1|2| — as well 
as from what is to be expected from LHC — , to 
lead to feasible precision tests of the Standard 
Model and trials of several New Physics scenarios, 
requires the knowledge of QCD matrix elements 
for their interpretation in terms of parameters of 
the Standard Model and its possible extensions. 
Unfortunately, the uncertainty on the theoreti- 
cal side in this interplay of experiment and the- 
ory in flavour physics predominantly originates 
from hardly computable long-distance effects of 
the strong interaction that confines quarks and 
gluons within hadrons. This potentially limits the 
impact of future experimental measurements on 
New Physics models and motivates calculations 
in lattice QCD, which is a powerful approach 
to reach a few-% theoretical error on those non- 
perturbative hadronic contributions. 

Still, some care is needed to obtain reliable re- 
sults for b-quark physics from a Monte Carlo eval- 
uation of the discretized Euclidean path integral. 
One has to keep under control simultaneously the 
finite-size effects and, particularly, the discretiza- 
tion effects, since the lattice spacing should not be 
larger than the Compton length of the b-quark. 
In practice, it is not possible to control both ef- 
fects by brute force numerical simulations such 
that dedicated methods have to be devised. 

While the numerical computations in lattice 
QCD necessarily involve approximations, one of 



the key features of the lattice approach is that all 
approximations can be systematically improved. 
For an overview of the different formulations of 
heavy quarks on the lattice that have been pro- 
posed in the literature and are being used to- 
day, and of results from the field of heavy flavour 
physics, which reflect some of these improvements 
by the small error bars quoted for many quanti- 
ties, I refer to the reviews of past Lattice Confer- 
ences |3I4I5I6I7] and references therein. 

1.1. Challenges 

Among the various considerable challenges one 
faces in an actual lattice QCD calculation on 
the theoretical and technical levels, let us only 
highlight the multi-scale problem, which is also 
particularly relevant in view of B-physics appli- 
cations. This is illustrated in Figure [T] There 
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Figure 1 . Large range of energy (/i ) scales in lat- 
tice QCD, where shaded areas refer to quark mass 
values (in the MS scheme) quoted by the Particle 
Data Group j^- -Red marks indicate the pion, the 
D- and the B-meson mass. 



are many disparate physical scales to be cov- 
ered simultaneously, ranging from the lightest 
hadron mass of « 140 MeV over tod ~ 2 GeV 
to tob ~ 5 GeV, plus the ultraviolet cutoff of 



1 



2 



Auv = of the lattice discretization that has 
to be large compared to all physical energy scales 
for the discretized theory to be an approximation 
to the continuum one. Moreover, the finiteness of 
the linear extent of space-time, L, in a numerical 
treatment entails an infrared cutoff Air = L^^ so 
that the following scale hierarchy is met: 

Air = L^^ ^ TO^, . . . , TTiD, me ^ — Auv ■ 

This implies L > 4/ m^r ~ 6 fm to suppress 
finite-size effects in the light quark sector and 
a < l/(2mD) ~ 0.05 fm to still properly resolve 
the propagation of a c-quark in the heavy sector. 
Lattices with L/a > 120 sites in each direction 
would thus be needed to satisfy these constraints, 
and since the scale of hadrons with b-quarks was 
not even included to arrive at this figure, it is ob- 
vious that the b-quark mass scale has to be sep- 
arated from the others in a theoretically sound 
way before simulating the theory. In Section [2] I 
briefly describe, how this is achieved by recours- 
ing to an effective theory for the b-quark. 

Another non-trivial task is the renormalization 
of QCD operators composed of quark and gluon 
fields, which appear in the effective weak Hamil- 
tonian, valid at energies far below the electroweak 
scale. Besides perturbation theory (see, e.g., 0), 
powerful non-perturbative approaches have been 
developed (and reviewed, e.g., in [10]), and I will 
come back to the non-perturbative subtraction of 
power-law divergences in the context of the effec- 
tive theory for the b-quark later. 

1.2. Perspectives 

As for the challenges with light quarks, it 
should only be noted that the condition L > 6 fm 
may be relaxed by simulating at unphysically 
large pion masses, combined with a subsequent 
extrapolation guided by chiral perturbation the- 
ory [11] and its lattice-specific refinements. 

Regarding the algorithmic side of a lattice QCD 
simulation, the Hybrid Monte Carlo [H] (HMC) 
as the first exact and still state-of-the-art algo- 
rithm has received considerable improvements by 
multiple time-scale integration schemes [13|14j . 
the Hasenbusch trick of mass-preconditioning [TSl 
1161 . supplemented by a sensible tuning of the al- 
gorithm's parameters |17| , and domain decompo- 



sition (DD) applied to QCD |18|19|20| . just to 
name a few. In addition, low-mode deflation |21j 
(together with chronological inverters has 
led to a substantial reduction of the critical slow- 
ing down with the quark mass in the DD-HMC. 

Finally, in parallel to the continuous increase 
of computer speed (at an exponential rate) over 
the last 25 years and the recent investments 
into high performance computing at many places 
of the world, the Coordinated Lattice Simula- 
tions [Ji (CLS) initiative is a community effort 
to bring together the human and computer re- 
sources of several teams in Europe interested in 
lattice QCD. The present goal are large- volume 
simulations with Ni — 2 dynamical quarks, us- 
ing the rather simple 0(a) improved Wilson ac- 
tion to profit from the above algorithmic devel- 
opments such as DD-HMC, and lattice spacings 
a = (0.08-0.05) fm, sizes L = (2-4) fm and pion 
masses down to m^r — 250 MeV, which altogether 
help to diminish systematic and statistical errors. 
Amongst others, the B-physics programme out- 
lined here is investigated on these lattices. 

2. Non-perturbative HQET 

Heavy Quark Effective Theory (HQET) at zero 
velocity on the lattice 24: offers a reliable solu- 
tion to the problem of dealing with the two dis- 
parate intrinsic scales encountered in heavy-light 
systems involving the b-quark, i.e., the lattice 
spacing a, which has to be much smaller than 
1/mb to allow for a fine enough resolution of the 
states in question, and the linear extent L of the 
lattice volume, which has to be large enough for 
finite-size effects to be under control (Figure [T]). 

Since the heavy quark mass (mb) is much larger 
than the other scales such as its 3-momentum or 
Aqcd 500 MeV, HQET relies upon a system- 
atic expansion of the QCD action and correlation 
functions in inverse powers of the heavy quark 
mass around the static limit (mt — > oo). The 
lattice HQET action Shqet at 0(l/TOb) reads: 

"''^J^x'^h {Do+6m- WkinD^ - WspinO-B} Tph , 

with i/'h satisfying P+tpi, = ^h, P+ = 

and the parameters Wkin and Wgpin being formally 

0(l/r7ib). At leading order (static limit), where 
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the heavy quark acts only as a static colour source 
and the light quarks are independent of the heavy 
quark's flavour and spin, the theory is expected to 
have ~ 10% precision, while this reduces to ^ 1% 
at 0(l/mb) representing the interactions due to 
the motion and the spin of the heavy quark. As 
crucial advantage (e.g., over NRQCD), HQET 
treats the 1/mb-corrections to the static theory 
as space-time insertions in correlations functions. 
For correlation functions of some multi- local fields 
O and up to 1/ rrib-corrections to the operator it- 
self (irrelevant when spectral quantities are con- 
sidered), this means 

{O) = (0),tat + {^kin(OOkin(a:))stat 

X 

+ Wspin(CCspin(2;))stat| , 

where (C')stat denotes the expectation value in 
the static approximation and Okin and Ospin are 
given by ■0h'-'^V'h and '0h''"BV'h- In this way, 
HQET at a given order is (power-counting) renor- 
malizable and its continuum limit well defined, 
once the mass counterterm Sm and the coeffi- 
cients Wkin and Wgpin are fixed non-perturbatively 
by a matching to QCD. 

Still, for lattice HQET and its numerical appli- 
cations to lead to precise results with controlled 
systematic errors in practice, two shortcomings 
had to be left behind first. 

1. ) The exponential growth of the noise-to- 
signal ratio in static-light correlators, which is 
overcome by a clever modification of the Eichten- 
Hill discretization of the static action [25] . 

2. ) As in HQET mixings among operators of 
different dimensions occur, the power-divergent 
additive mass renormalization Sm ^ ja already 
affects its leading order. Unless HQET is renor- 
malized non-perturbatively [26], this divergence 
— and those ^ g^jo^ arising at 0(l/mb) — im- 
ply that the continuum limit does not exist owing 
to a remainder, which, at any finite perturbative 
order |27l28j . diverges as a — 0. A general so- 
lution to this theoretically serious problem was 
worked out and implemented for a determination 
of the b-quark's mass in the static and quenched 
approximations as a test case [22] . It is based on a 
non-perturbative matching of HQET and QCD in 



finite volume. Applications of this strategy to the 
determination of the b-quark mass and (a subset 
of aU) HQET parameters at 0(l/mb) |30I31| . to 
a study of the Bg-meson spectrum !32J and to a 
computation of the Bg-meson decay constant }33) 
were realized in the quenched approximation by 
our collaboration and have been extended to the 
more realistic Nf — 2 situation [34 35 3 6 ,37] . 

3. The b-quark mass via HQET at 0(l/mb) 

We first note [3H| that in order not to spoil the 
asymptotic convergence of the series, the match- 
ing must be done non-perturbatively — at least 
for the leading, static piece — as soon as the 
l/mb-corrections are included, since as mb — > oo 
the perturbative truncation error from the match- 
ing coefficient of the static term becomes much 
larger than the power corrections ~ Aqcd / m-b of 
the HQET expansion. 

QCD ! HQET 



Li I Li L2 L2 L- 




Si I 52 S3 Si Ss 



Figure 2. Idea of lattice HQET computations 
via a non-perturbative determination of HQET 
parameters from small-volume QCD simulations. 
For each fixed Li, the steps are repeated at 
smaller a to reach the continuum limit. 

In the framework introduced in 29i, match- 
ing and renormalization are performed simultane- 
ously and non-perturbatively. Let us here explain 
the general strategy, illustrated in Figure [2] for 
the sample application of calculating the b-quark 
mass. 5*1: Starting from a finite volume with 
Li w 0.5 fm, one chooses lattice spacings a suffi- 
ciently smaller than l/wb such that the b-quark 
propagates correctly up to controllable discretiza- 
tion errors of order a^. The relation between the 
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renormalization group invariant (RGI) and the 
bare mass in QCD being known, suitable finite- 
volume observables $fc(Li, Af^) can be calculated 
as a function of the RGI heavy quark mass, Mh, 
and extrapolated to the continuum limit. 6*2: 
Next, the power-divergent subtractions are per- 
formed non-perturbatively by a set of matching 
conditions, in which the results obtained for $fc 
are equated to their representation in HQET. 
At the same physical value of Li but for reso- 
lutions Li/a — 0(10), the previously computed 
heavy-quark mass dependence of $fc(Li,Mh) in 
finite- volume QCD may be exploited to determine 
the bare parameters of HQET for a w (0.025 — 
0.05) fm. 53: To evolve the HQET observables 
to large volumes, where contact with some phys- 
ical input from experiment can be made, one 
also computes them at these lattice spacings in a 
larger volume, L2 = 2Li. The resulting relation 
between $fc(Li) and $fc(L2) is encoded in associ- 
ated step scaling functions denoted 
By using the knowledge of ^k{L2, M^) one fixes 
the bare parameters of the effective theory for 
a ~ (0.05 — 0.1) fm so that a connection to lattice 
spacings is established, where large-volume ob- 
servables, such as the B-meson mass or decay con- 
stant, can be calculated. This sequence of steps 
yields an expression of tob, the physical input, as 
a function of Mh via the quark mass dependence 
of <I>fe(Li, Mh), which eventually is inverted to ar- 
rive at the desired value of the RGI b-mass within 
HQET. The whole construction is such that the 
continuum limit can be taken for all pieces. 

3.1. Computation of HQET parameters 

Following the strategy sketched above and ap- 
plied to the quenched case in [31 , the determina- 
tion of the parameters of the HQET Lagrangian 
and of the time component of the isovector axial 
current is performed within the Schrodinger func- 
tional, i.e., QCD with Dirichlet boundary condi- 
tions in time and periodic ones in space, where 
suitable matching observables such as finite- 
volume meson energies and matrix elements, can 
be readily defined. Relativistic quarks are sim- 
ulated as clover-improved Wilson fermions with 
Nf — 2 dynamical quarks; for the static quark we 
use the so-called HYP 1/2 actions pS] . 



We introduce observables ^k=i,...,5 casted into 
a vector $ = <i)Q'^'^ , where in the continuum and 
large volume limits, the first two are proportional 
to the meson mass and to the logarithm of the 
decay constant, respectively, while ^3 is used to 
fix the counterterm of the axial current and $4.5 
for the determination of the kinetic and magnetic 
terms in iShqet- The continuum extrapolations 
of $1,2 in the small QCD volume {Li « 0.5 fm, 
Si in Figure [2| , for nine values Mh = M of the 
RGI heavy quark mass from the charm to beyond 
the bottom region [36j, are shown in Figure [sj 
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Figure 3. Continuum extrapolation of the Rnite- 
volume observables <&! and $2, where for $1 we 
have included the error (cross on the left) stem- 
ming from the renormalization of the quark mass. 

When the effective theory is simulated in the 
same physical volume (5*2 in Figurej2]), a set of 
matching conditions for lima_s.o (ii, M, a), 

<l>'^^^{Li,M) - 77(^1, a) + ^iLi,aMM, a) , 

is imposed; the r.h.s. represents the heavy quark 
mass expansion of the ^Q*^^ at 0(l/TOb)- Hav- 
ing computed 77 and (j) from these simulations for 
different values of a, the matching equations de- 
termine the set of parameters a)(M, a); e.g., in the 
simple case of the static meson mass, and up to 
a kinematic constant, rj is the static energy, (j) a 
constant and uj the bare static quark mass. After 
step scaling to L2 = 2Li, the observables in this 
volume are now obtained, thanks to the parame- 
ters uj{M, a) fixed by the previous step, as 

$(L2,M,0) = lim [r]{L2,a) + <j){L2,a)Cj{M,a)] , 

a— 

and the continuum limit can be taken, since the 
power divergences in HQET cancel out here. Fig- 
ure |4] depicts examples of corresponding contin- 
uum extrapolations in the static approximation, 
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and the results for observables sensitive to the 
1 /mb-corrections are of similar quality [37) . 
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Figure 4. Continuum extrapolation of the static 
approximation of <f>i and $2 in the volume of ex- 
tent L2- Red (blue) symbols refer to the HYPl 
(HYP2) discretization of the static propagator. 

Finally, the HQET parameters to be employed 
in the large volume, Lqoi Siie estimated from 6*3: 

oj{M,a) ^ (b-\L2,a)mL2,M,0) - r^iL2,a)] . 

3.2. Preliminary large-volume results 

To apply the non-perturbative matching results 
to calculate the b-quark mass, we write down the 
HQET expansion (to first order in l/mb) of me 
in terms of HQET parameters and energies as 



me = mba 



C^ki„i?^'"+Wspi„i?^P'". (1) 



HQET energies and matrix elements have been 
extracted from measurements on a subset of con- 
figuration ensembles produced within CLS [23] 
solving the Generalized Eigenvalue Problem [35] , 
which allows for a clean quantification of sys- 
tematic errors from excited state contaminations. 
So far, only a single lattice spacing a « 0.07 fm 
(/3 = 5.3) has been analyzed so that the size of 
discretization eS^ects can not be assessed yet. Fig- 
ure [5] shows elements of the computations in Loo- 
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Figure 5. Left: Comparison of plateaux of static 
energies at (3 = 5.3 to earlier quenched results. 
Right: Graphical solution of ([I]) in static approx- 
imation; M = Mh is the RGI heavy quark mass. 



Upon chiral extrapolation in the light quark 
mass and including a conservative uncertainty in 
the lattice scale (ro = 0.475(25) fm [iD]), we quote 
as our preliminary result for the b-quark's mass 
in HQET at 0(l/mb) for the A^f = 2 theory: 

mb^(mb) = 4.276(25),„(50),tat+rcnorm(?)aGeV. 

The first error states the scale uncertainty, while 
the second covers the statistical errors of HQET 
energies, the chiral extrapolation uncertainty and 
the error on the quark mass renormalization en- 
tering the small- volume QCD part of the compu- 
tation (Si). More details are found in [37] . 

For comparison, we cite the previous iVf = 
HQET result m^(mb) = 4.320(40)^„ (48) GeV 
by our collaboration [30] and the recent sum-rule 
(Mb) = 4.163(16) GeV [H]. 



determination, 



4. Outlook 

The non-perturbative treatment of HQET in- 
cluding l/mb-terms can lead to results with un- 
precedented precision for B-physics on the lattice. 
It also greatly improves our confidence in the use 
of the effective theory. Our project to extract 
from N{ = 2 lattice simulations relevant quanti- 
ties for B-phenomenology within HQET is well 
advanced. While the non-perturbative matching 
of HQET with QCD through small-volume sim- 
ulations is almost done, the evaluation of HQET 
energies and matrix elements has started recently 
on the CLS ensembles, but still awaits a bet- 
ter control of the cutoff effects. Our first results 
for the b-quark mass in HQET at 0(l/mb) are 
promising, and further applications of the once 
determined HQET parameters to calculate the 
B-meson decay constant, the spectrum of heavy- 
light mesons and the form factors of the B — >■ tt 
semi-leptonic decay are expected in the future. 
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